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^ ■ Abstract 

^ ■ The two-nucleon spectral function in nuclear matter is studied using Correlated 

'sj" , Basis Function perturbation theory, including central and tensor correlations pro- 

\ duceded by a realistic hamiltonian. The factorization property of the two-nucleon 

Q"^ ' momentum distribution into the product of the two single nucleon distributions 

. shows up in an analogous property of the spectral function. The correlated model 

yields a two-hole contribution quenched whith respect to Fermi gas model, while the 
peaks acquire a quasiparticle width that vanishes as the two momenta approach kp- 
> . In addition, three-hole one-particle and more complicated intermediate states give 

\ rise to a background, spread out in energy and absent in the uncorrelated models. 

^ ' The possible connections with one- and two-nucleon emission processes are briefly 

^ ■ discussed. 

> : 
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1 Introduction 



The two-nucleon spectral function, yielding the probability to remove two 
nucleons with momenta ki and k2 from nuclear matter leaving the residual 
system with excitation energy E, is defined as (see, e.g., ref.[l]): 

P(ki, k2, E) = Y. mai^ai^\n{A - 2))\^5{E - E^ - 2) + Eo{A)) , (1) 

n 

where |0) is the A-particle nuclear matter ground state with energy Eq{A), 
\n{A — 2)) denotes a (A — 2)-nucleon intermediate state with energy En{A — 2) 
and ajj is the usual creation operator. 
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The two-micleon spectral function carries relevant information on the short 
range structTirc of the nuclear medium and nucleon-nucleon (NN) correlations. 
Therefore, a considerable effort aimed at extracting experimental information 
on P(ki, k2, E) is currently being undertaken. In this context, a very impor- 
tant role is played by electron-nucleus scattering experiments (a number of 
theoretical and experimental topics in the field of electron-nucleus scattering 
are reviewed in ref.[2]). 

Unambiguous evidence of strong correlation effects has been provided by single 
nucleon knock out (e, e'p) reactions, showing that in a nucleus only about 70 % 
of the nucleons are in states of low momentum and low removal energy, that 
can be described by a mean field theory. The remaining 30 % of the nucleons 
belong to strongly correlated pairs, whose occurrence is mainly to be ascribed 
to the one-pion-exchange tensor force and to the repulsive core of the NN 
interaction. 

The available (e, e'p) data at low missing momentum and low missing en- 
ergy clearly show the depletion of the occupation probabilities of single par- 
ticle states predicted by the nuclear shell model (for a recent review see, e.g., 
ref.[3]), thus providing a somewhat indirect measurement of correlation effects. 
Complementary information will soon be available from single nucleon knock 
out experiments specifically designed to investigate the kinematical region 
corresponding to large missing momentum and large missing energy, where 
correlation effects are believed to be dominant [4]. 

Over the past decade, the availability of a new generation of high energy 100 % 
duty-cycle electron beams has made it possible to carry out double coincidence 
(e, e'NN) experiments, in which two knocked out nucleons are detected. In 
principle, these experiments may allow for direct measurements of correlation 
observables, such as the momentum space two-nucleon distribution or the two- 
nucleon spectral function. However, as pointed out in ref. [5], extracting the 
relevant dynamical information from the measured (e, e'NN) cross section 
may turn out to be a challenging task, requiring both a careful choice of the 
kinematical setup and a quantitative theoretical understanding of the final 
state interactions (FSI) of the knocked out nucleons. 

Regardless of the difficulties involved in the interpretation of the experimental 
data, reliable theoretical calculations of the two-nucleon emission cross section 
carried out within the Plane Wave Impulse Approximation (PWIA), in which 
all the nuclear structure information is contained in the two-nucleon spectral 
function, have to be regarded as a minimal starting point, which can provide 
guidance for the optimization of future experiments. 

The basic assumptions underlying the PWIA scheme are that i) the exchanged 
virtual photon couples to either of the two outgoing nucleons and ii) FSI effects 
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are negligible. The PWIA cross section of the process in which an electron of 
initial energy Ei is scattered into the solid angle fig with energy Ej = Ei — u, 
while two nucleons of kinetic energies Tp and are ejected into the solid 
angles flp and fl'^, respectively, takes the simple factorized form 

du^dnjnpdTpdn'^dT; = + ^^^'^^^ + m)a,^F,AP^, E^) , (2) 



where Uem describes the structure of the electromagnetic vertex, m is the 
nucleon mass and p and p' denote the momenta of the detected nucleons. The 
missing momentum p^ and missing energy E^ are defined as 

P.. = q-P-P', (3) 



q being the momentum transfer, and 

E^ = uj-Tp-T'-Tr, (4) 



where Tr is the kinetic energy of the recoiling [A — 2)-particle system. The 
function Fpp'(pm, Em) appearing in eq.(2) can be written in terms of the two- 
nucleon spectral functions P(p — q, p', and P(p, p' — q, Em)- 

Pioneering studies of the PWIA two-nuclcon emission cross section of the 
^^C(e, e'pp) reaction are discussed in ref.[6], whereas a G-matrix perturba- 
tion theory calculation of the ^^O two-proton spectral function is described in 
ref.[7]. 

In this paper we discuss a calculation of the nuclear matter P(ki, k2, E) per- 
formed using Correlated Basis Function (CBF) perturbation theory. Our the- 
oretical approach, in which the effects of the nonperturbative components of 
the NN interaction are incorporated in the basis function using a variational 
approach, has proved to be particularly suited to describe quantities which 
are strongly affected by NN correlations, such as the single nucleon spectral 
function P(k, E)[8], the nucleon momentum distribution n(k)[9] and the off- 
diagonal density matrix p(ri, r']^)[10]. 

Following ref.[8], we evaluate the dominant contributions to P(ki,k2,£') at 
the zero-th order of CBF, as well as the perturbative corrections that have 
been shown to be relevant in the calculation of P(k, E). 

CBF theory of infinite nuclear matter is built on the set of correlated states: 

\N)^S X{F{t,3) \N)fg. (5) 
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obtained by applying a symmetrized product of two-body correlation opera- 
tors, F{i,i), to the Fermi gas states \N)fg. The structure of F{i,j) is similar 
to that of the NN interaction: 

F{i,j)^j:f-\n,)0^-\i,j). (6) 

n 

The operators 0^'^\i.j) include four central components (1, (ai-aj), (Ti-Tj), [oi- 
(Tj){Ti ■ Tj)) for n = 1,4 and the isoscalar and isovector tensor ones Sij and 
Sij{Ti ■ Tj) for n = 5,6. Additional spin-orbit components are sometimes in- 
troduced, but we will neglect them in this work. With transparent and widely 
employed notation, the components are also denoted as c {n — 1), a, t and t 
(tensor) . 

Realistic correlated basis states are expected to be close to the eigenstates of 
the hamiltonian. An efficient recipe to choose the correlation functions consists 
in their variational determination by minimizing the expectation value of the 
hamiltonian in the correlated ground state: 

K = m\o) . (7) 

The scalar (or Jastrow) component /("=i)(r) heals to unity for r — > oo, whereas 
/("^^^(r) — >■ 0. The tensor components have the longest range, since they are 
related to the one-pion exchange potential and F\i, j)HijF[i, j) {Hij is the 
hamiltonian associated with a two-nucleon pair) is a well-behaved operator 
healing to Hij at large interparticle distances, like the G-matrix. The results 
presented in this work have been obtained using the Urbana V14 potential 
supplemented by the TNI three-nucleon interaction [11]. 

At the zero-th (or variational) order of CBF, P^arO^i, ^2, E) is obtained using 
the correlated states (5) in eq.(l). The variational estimate is then corrected by 
inserting perturbative corrections in the correlated basis. This procedure has 
been already adopted for many nuclear matter properties and it has been par- 
ticularly successful in describing inclusive electromagnetic responses at both 
intermediate and high momentum transfers (see, e.g., ref.[2]). 

In Section 2 we will discuss the relationships between the two-nucleon spectral 
function, the two-nucleon momentum distribution, and the two-body density 
matrix, p(ri, r2; r^, r2). Section 3 will present theory and results for the vari- 
ational two-nucleon spectral function, as well as the main CBF perturbative 
corrections, evaluated with realistic interactions. Finally, Section 4 is devoted 
to conclusions and perspectives. 
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2 Spectral functions, momentum distributions and density matri- 
ces 

The nuclear matter two-body density matrix p(ri, r2, r'^, r'^ is defined as 

p(ri,r2;r;,r'2)=A(A- 1) 

/ rfVg ■ ■ ■ d^TA ra, ■ ■ ■ , r^)^o(r;, v'^, . . . , r^) 

/#ri...#^^|*o(ri,r2,...,r^)|2 ' ^> 

^o(ri, r2, . . . , r^) being the ground state wave function. For ri = r'^ and r2 = 
r'21 the density matrix reduces to the two-nucleon density distribution p(ri, r2), 
yelding the joint probability of finding two nucleons at positions ri and r2 in 
the nuclear matter ground state. 

The two-nucleon momentum distribution, i.e. the probability that two nucle- 
ons carry momenta ki and k2, is related to the two-nucleon spectral function 
P(ki, k2, -B), defined in eq.(l), through 

n(ki,k2) = j dEP{ki,k2,E) = (0|aj,^aj,^ak,aki|0) . (9) 

From Eqs.(8) and (9) it follows that n(ki, k2) can be written in terms of the 
density matrix according to 

n(ki,k2) = J rfVidV2dMt^Me'''"'""'e"'"'"^^'p(ri,r2;r;,r'2) , (10) 

where u denotes the degeneracy of the system (in symmetric nuclear matter 
v = 4), Q is the normalization volume, rw = ri — ri/ and r22' = r2 — r2'. 

It is well known [12] that in a Fermi liquid with no long range order the two- 
nucleon momentum distribution factorizes according to 

n(ki,k2) = n(ki)n(k2) + oQ) , (11) 

where n(k) is the single nucleon momentum distribution, defined as 

n(k) = (0|aj,ak|0) . (12) 

The n(k) resulting from the nuclear matter calculation of ref.[9] is shown in fig. 
1. It exhibits a discontinuity at |k| = kp (the Fermi momentum kp is related 
to the density p — A/Q through kp = (Gtt^p/z/)^/^), and a tail, extending to 
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very large momenta, reflecting the structure of the nuclear wave function at 
small interparticlc distance. Eq.(ll) implies that in the A ^ oo limit n(ki, 
contains the same dynamical information carried by n(k). 

Some insight on the role of short range correlations can be obtained rewriting 
n(ki,k2) in terms of the relative and center of mass momenta of the pair, 
defined as q = (ki — k2)/2 and Q = ki + k2, respectively, and studying the 
quantity 

n^ez(q) = 47r|qp /d^Q n 



Q 



+ q 



n 



Q 



- q 



(13) 



Fig. 2 illustrates the behaviour of n^e^q) saturation density (corresponding 
to kp — 1.33 fm""^), evaluated using the momentum distribution of ref.[9] 
(solid line) , compared to the prediction of the Fermi gas model (dashed line) . 

The inclusion of NN interactions leads to a quenching of the peak, located 
at |q| ~ 0.7-0.8 fm~^, and to the appearance of a sizeable tail at |q| > kp- 
It is interesting to single out the contributions to nrei(q) corresponding to 
ki,k2 < kp, ki < kp and k2 > kp, and ki,k2 > kp. They are represented 
in fig. 2 by diamonds, squares and crosses, respectively (note that the results 
shown by the crosses have been enhanced by a factor 10). The ki,k2 < kp 
component provides about 72 % of the normalization of nrei(q) and becomes 
vanishingly small at |q| > kp. On the other hand, the contributions coming 
from pairs in which at least one of the nucleon momenta is larger that kp 
are much smaller in size (the curves marked with squares and crosses yield 
26 % and 2 % of the normalization, respectively) but extend up to |q| ~ 4 
fm^^. Similar results have been recently obtained [13] using the momentum 
distribution resulting from a G-matrix perturbation theory calculation [14] . 



3 CBF calculation of the two-nucleon spectral function 



The propagation of two hole states in nuclear matter is described by the hole- 
hole part of the two-body Green's function[l], 

Ghh{^i,k2,E) = {0\al^al- ^ ^ ak^QkilO) , (14) 



where H is the nuclear matter hamiltonian. Ghh is defined for e(= —E) < 2ep, 
ep being the Fermi energy. 
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The two-nucleon spectral function, given in eq.(l), is straightforwardly related 
to the imaginary part of Ghh by 

P{k,,k2,E) = -^Ghh{kiM,E) . (15) 



CBF perturbation theory employes the set of correlated states given in eq.(5). 
Note that the basis states are non orthogonal to each other. Orthogonahza- 
tion can been implemented using the procedure proposed in ref.[15], which 
preserves the diagonal matrix elements of the hamiltonian between correlated 
states. This technique has been succesfuUy applied to the calculation of the 
one-body spectral function in ref.[8] (hereafter denoted I). 

The perturbative expansion is obtained separating the nuclear hamiltonian H 
into two parts: H = Hq + Hj, and using standard Rayleigh-Schrodinger type 
expansions. The unperturbed and interaction terms, denoted by Hq and Hj, 
respectively, are defined through 

{N\Ho\M) = 5nm{N\Ho\N) = SnmE"^ , (16) 

{N\Hi\M) = (1 - 6nm){N\Hj\M) . (17) 

The expectation value in eq.(14) is calculated by expanding 

^ Eni-r [{H o - E^,r' {Hi - AEo)Y |0) 



|(0|0)|V2 



where AEq — Eq — Eq. Perturbative corrections to the intermediate states 
|A^) arc taken into account carrying out an analogous expansion in powers of 
{Hj — AEq) for the propagator appearing in eq.(14), as discussed in I: 

1 _ 1 

{H - Eo - E - ir])~ {H - El - E - ir]) 

xi:(-r{w-A^o)(,,3^^}".(i9, 

We will now describe the structure of the n — (i.e. variational) and higher 
order terms of the CBF perturbative expansion of P(ki, k2, £^). 
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3. 1 The variational spectral function 



At the lowest order in the perturbative expansion the spectral function is given 

by 

P,„.(ki, k2, E) = E l(0|4«lcJr^(^ - 2))\'6{E - E^^A - 2) + El{A)) , (20) 

n 

where the correlated states (5), and their energies Eq{A) and E^{A — 2) are 
used. 

The main contributions to the sum in eq.(20) come from correlated two- hole 
(2/i) and three-hole one-particle (Ship) intermediate states. In the uncorre- 
lated Fermi gas only 2h states contribute to P(ki,k2, E), leading to the well 
known result 

P2hM^iM,E)^e{kF - \k,\)e{kF - |k2|) 

where Q{x) denotes the usual step function. The presence of correlations, 
besides changing the structure of the 2h contribution in the way we will discuss 
below, allows for nonvanishing contributions from other intermediate states, 
with a consequent quenching of the 2h peak since part of the strength is moved 
to higher excitation energies. 

The variational 2h spectral function is given by: 



X Siel + el + E) , (22) 

where e)^ = (h|i/|h) — E^ is the variational single particle energy and the 2h 
overlap matrix element, ^k^ kj) 

<'£ = (0|aLaL|hi,h2). (23) 



The correlated Ih overlap matrix element, — (0|^k|h), was computed in 
I using cluster expansion techniques. The same method is used for the 2h 
overlap, with the result that only unlinked cluster diagrams {i.e. diagrams 
where the points reached by the /ci^2-lines are not connected to each other by 
any dynamical, P — 1 and / — 1, or statistical correlations) contribute. The 
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linked diagrams contribution turn out to be of order 1/A, and therefore vanish 
in nfinite nuclear matter. Note that the factorization property of the two- 
nucleon momentum distribution, illustrated by eq.(ll), is also a consequence 
of the fact that linked diagrams do not contribute in the A— > oo limit. The 
definitions, as well as the technicalities entering the cluster expansion method 
can be found in I and in the related references and will not be repeated here. 
Because of the factorization property, in nuclear matter we obtain 

<;£ = <^>ki<^5h,-k.<^h.-k. , (24) 



the CBF expression of being given in I (eq.(A.19)). 

The correlated 2h spectral function retains the Fermi gas delta shaped peak. 
However, correlations quench the peak itself (via the 2h overlap matrix ele- 
ment) and move it to = — e)^^ — ej^^. 

The Ship contribution reads 



hl,h2,h3,pi 

X S{el + el + el - e^ + E) , (25) 
the overlap matrix element being given by 

^h.,h„h3,p, ^ (olaj^^aLjhi, h2, h3, pi) . (26) 

In the limit A— >^ oo the contributions of fully linked diagrams vanish in the 
Ship overlap as well. As a consequence, 

^hi,h2,h3,pi _ PI /r,-7X 

In the other cases we obtain 

^hi,h2,h3,pi _ \^ ds'^l<T>'^2,h3,pi r /qo\ 

*ki<fep-,k2>feF - Z^*ki*k2 <)hi-ki , i^ej 



and 



<<t':k2"<.. = E [<^<"'""^^h,-k, + - k^] , (29) 



where the sums include all permutations of the hj indices, while the structure 
of the 2hlp overlap ^j^i''^^,?! _ ^Qlaj^-lhi, h2, pi) is discussed in I (Appendix 
B). 
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Unlike the 2h contribution, P3hip,v is no longer delta shaped, but it is rather 
spread out in energy. It may be interpreted as a background contribution to 
be added to the two-quasiparticle part of the spectral function. 

The next intermediate state to be considered is the 4h2p one, whose contribu- 
tion may be written in terms of products of two squared 2hlp overlap matrix 
elements. However, it is numerically negligible with respect to the Ship term 
(several order of magnitudes smaller) and can be safely neglected. 

The above intermediate states completely exhaust the momentum distribution 
sum rule in nuclear matter at the variational level: 



^ j dEP^ar{^l,\i2,E) . (30) 



3.2 Perturbative corrections 



The admixture of m-hole n-particle correlated states in \N) originate CBF 
perturbative corrections. As far as the ground state is concerned, they corre- 
spond to the n>l terms in eq.(18). 

Following the strategy developped in 1, the perturbative corrections are clas- 
sified as dPcR{^i,i^2, E) and 6PjNT(ki,^2, E), coming from ground and in- 
termediate state corrections, respectively. The total spectral function is then 
given by 

P{k,M:E)=P,ar{^lM:E) + 6PGR{k,M:E) + dPiNT{^,M:E) .(31) 



Here we consider 2h2p admixtures to the ground state and 2hlp admixtures 
to the Ih intermediate states. The first ones contribute for all |k| values (both 
below and above kp), and their contribution can be expressed in terms of the 
correlation self energy E'-^'^(k, E), whose imaginary part is given by 

Qj:^0{k,E) ^lEmH\kp2hM\'5{E + e;^-el-el) , (32) 

at |k| > kp and E < and 

QE^^iKE) ^lj:mH\PiP2^i^2)mE + el-el-el) , (33) 
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at |k| < kp and E > e^p. The self energy is computed using correlated states 
and retaining two- and three-body separable contributions in the cluster ex- 
pansion of the relevant matrix element [17]. 

The 2hlp intermediate state contribution involves the polarization self energy 
E'^'^(k, E). The corresponding imaginary part reads 

3S^«(k, = I j:\{k\H\p,h,h,)\'6{E + e;^ - el - elj , (34) 

at |k| < kp and E < Cp and 

»E^O(k,E) = lY:\{^\H\PiP2h2)mE + el-el-e;j , (35) 



at |k| > kp and E > Cp. 

The effects of these corrections on the spectral functions at |ki|, |k2| < kp are: 
i) a quenching of the variational two-quasiparticle strength and a shift of its 
position by a quantity Se^^ + Se^^{k), with 

5,co_l 7..^S^^(k,^) 



from ground state corrections, and ii) & broadening of the delta-like peak 
itself, roughly proportional to '^T.^'-^iki, E — — e^J -|- 'ts'E^'^{]s.2, E — —e^^), 
from intermediate state admixtures. 

The expression for the total 5P at |ki|, |k2| < /cf reads 



K;£i'+sr(ki,i?)+sr^(k2,s) 



71 ^ 

X [a{ki, k2){-el^ - e^^ - E) - Sck^ - 5ek^ 

- ^ (E^^(ki, E) - E^°(ki, E = -el) + 1^2) 

-t (S^^(ki,i?) + 1^2) 



(37) 



where a{ki, k2) = 1 — Sn2{ki) — Sn2{k2) — Sn2{ki) — Sn'2{k2). The quantities 
5n2{k), 5n2{k) and Ef'^(k, £■) are defined in I. 



Fig. 3 gives two examples of the two-nucleon spectral function below kp. The 
figure shows the variational one-hole {Pih,v) and two-hole {P2h,v) peaks and the 
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total spectral function (solid line). The one-nucleon spectral function below 
kp can be separated into a single particle and a correlated part [18]. The 
former originates from Ih intermediate states and their admixtures, the latter 
from n-hole {n — l)-particle states, and it would be strictly zero in absence 
of correlations. The corresponding contributions for the two-nucleon spectral 
functions arc shown in the figure as dashed and dot-dash lines, respectively. 
The correlated part has been enhanced by a factor 10^, while the single particle 
and total ones by a factor 10. The figure also gives the barely visible variational 
correlated part (dotted lines). 

Ground state 2h2p corrections give the leading contributions to SP when one 
of the momenta is above kp. So, for |ki| < kp and |k2| > kp, we obtain 



Spk^E+^l (38) 
where ^(k, E) is again given in I. 

The spectral functions for two such cases are shown in fig. 4, together with the 
one-hole variational peaks corresponding to the momenta below kp. The whole 
spectral function now comes from the correlated part. The figure gives the 
total (sohd lines) and variational (dotted lines) P, as well as the perturbative 
corrections (dashed lines). 

Corrections at |ki|, |k2| > kp are quadratic in the self energies and have not 
been considered. 



4 Summciry and conclusions 



We have carried out a calculation of the nuclear matter two-nucleon spec- 
tral function P(ki,k2,-E') within the framework of CBF perturbation theory. 
The zero-th order approximation, corresponding to a variational estimate, has 
been supplemented with higher order corrections associated with both 2h2p 
admixture to the ground state and 2hlp admixtures to the Ih states. 

The results show that the inclusion of NN interactions produces drastic changes 
in the behavior of P(ki,k2,-E), with respect to the predictions of the Fermi 
gas model. The peaks correponding to the single particle states get quenched 
already at zero-th order, and aquire a width when higher order corrections are 
included. Due to strong short range NN correlations, the strength removed 
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from the quasi particle peak is pushed to large values of momenta and re- 
moval energy, giving rise to a broad background. 

The factorization property exhibited by the two-nucleon momentum distri- 
bution of a normal Fermi liquid has been found to hold at the level of the 
amphtudes entering the calculation of P(ki, k2, E) as well. As a result, all the 
linked cluster contributions to the two-nucleon spectral function vanish in the 
infinite nuclear matter limit. 

In principle, the results discussed in this paper may be used to estimate cor- 
relation observables, related to the small components of P(ki,k2,-E), whose 
measurement represents the ultimate goal of the (e, e'NN) experimental pro- 
grams. However, it has to be emphasized that, in order to make fully quanti- 
tative predictions, our approach should be extended to finite systems, working 
along the line proposed in refs.[19,10]. The analysis of the A dependence of 
P(ki,k2,i?) would also provide very interesting information on the relevance 
of the 1/A linked contributions and the applicability of the factorized approx- 
imations for the two-nucleon spectral function and momentum distribution. 
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Fig. 1. Single nucleoli momentum distribution in nuclear matter, evaluated in ref.[9] 
using CBF perturbation theory. 
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Fig. 2. The distribution nrei{<l) defined in eq.(13). The solid and dashed fines 
represent the results obtained using CBF perturbation theory and the Fermi gas 
model, respectively. The meaning of the other curves is explained in the text. 
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